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MODELS AND MAPS OF THE MAIN FIELD 
by 

C. GIRE and J. DUCRUi:: 


OBJECTIVE : 

The first objective of our investigation is to obtain some models of the 
main field. The difference between these models comes from the difference in 
the procedures of data selection and depends on the aim the model is intended 
for : regional studies necessitates locally higher density of data^for downward 
continuation the reduction of external variations must be as careful as possible. 
We present here refinements of our preliminary models MAGPl and MAGP2 (see 
our progress reports entitled "Preliminary models of the core field") but 
these models are not yet definitive though they compare fairly well with those 
obtained by other teams (see for example Lange 1 , 1981). 

BACKGROUNDS : 

Our models MAGPl and MAGP2 were found to have drawbacks coming from 
insufficient reduction of the external variations and from the inhomogeneity 
in the distribution of the measurement points. New analysis has been performed 
to overcome these shortcomings. The quality of the new models had to be charac- 
terized by the average r.m.s. of the differences measurement-model as the 
facility of having profiles and maps rapidly drawn is not yet available to us. 

This generation of models has the same general characteristics as the 
preceding one. The raw data have been averaged over 2.5s (this is approximately 
18 km along the satellite track). To shorten the computations only one point 
out of sixty has been retained and the points above 60° of latitude have been 
discarded to lower the influence of the external perturbations. The component 
weighting has been kept to the same value : Px “ 1.2, Py ■ 0.6 and p^ ■ 


0 . 8 . 


We also keep the value N » 1 1 for the maximum order of the spherical harmonic 
development, this value being probably overestimated with respect to the 
precision obtained for the higher harmonic coefficients. 

RECENT ACCOMPLISHMENTS : 

We received the CHROFIN tapes from November 2”*^ 1979 to May 19^^ 1980. 

Four operations have to be made before having data ready for use in the spherical 
harmonic analysis. First translation from IBM binary code to CDC usable code, 
second averaging the measurements over a 2.5 second internal, third interpola- 
ting the satellite position and fourth determining the geomagnetic external 
activity index. This processing has been applied to the tapes containing the 
data from December 2"^*^ 1979 to April 4*-*^ 1980 (CHROFIN tape number 8 to 22). 

In five files problems occur making the whole file unreadable (see arrows 
in Table 1). We have not found yet a solution to these problems. 

spherical harmonic analysis : 

We tested the influence of the polar regions on the quality of the models. 

We choosed a small subset of the data (from 12/02/79 to 12/06/79) containing 

752 points having a magnetic activity index equal to 0 or 1 (for the meaning 

of this index see our progress report entitled "Separation of internal and 

external fields : a new technique of data screening"). Data points with a 

latitude larger than 60“ were weighted with a weight Pq varying between 0 and 

0.1. The corresponding average r.m.s of the residues are displayed in Table 2. 

These results show no dramatic influence of the polar regions so we decided to 

discard these regions for our further analysis. However it has to be noted 

the*’ *’he g*** and h™ coefficients are very stable for n< 7 and seems to be Po- 
n n 

dependant for n>7. 

We prepared three models using the data from 12/02/79 to 01/10/80. The 
first one called MGSTOl used the same time interval of data as MAGPl (from 


12/01/79 to 12/21/79) with 2436 data points selected according to their index 

(0 or 1). The second one called MGST02 used no selection according to the e ter- 

nal magnetic index and used 7247 measurements. The third one, MGST03, used 

1996 data points with an index 0 or 1 belonging to the period 12/22/79 to 

01/10/80. The average r.m.s of the residues of these models are displayed in 

Table 3. We decided to retain model MGST03 because the corresponding data points 

distribution is more uniform than for the MGSTOl model. The g'^ md h™ coefficient 

tl n 

of the MGST03 model are shown in Table 4. The corresponding maps for the ele- 
ments F,D,I and Z have been drawn (see figures 1 to 4). 

anomaly maps : 

Tentative anomaly maps have been computed but are not yet satisfactory. 

We used S^xS** averaging which is too coarse and the external variation reduction 
needs to be refined. Further attempts are to be made (see next section). 

secular variation : 

The secular variation is smaller than 0.5nT/day all over the surface of 
the Earth. Then we decided to neglect this variation as a first approximation 
the resulting error being lower than 10 nT for our MGST03 model. We undertook 
a detailed study of the secular variation (Ducruix et al, 1980; Le Mouel et al, 
1982) and we will apply our knowledge to the preparation of a precise spherical 
harmonic model of the secular variation (see next section). 


FUTURE EMPHASIS : 

Our next step is the obtaining of precise anomaly maps. The critical point 
for this is the reduction of the external variation which is starting now using 
manual screening. We are also developping a algorithm of automatic screening 
based on pattern recognition but this will probably take some time to be 


achieved. 
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The secular variation will be described and taken into account using our 
model and a new method which will be displayed in our final report. We begun 
the analysis of the motion at the core surface corresponding to the secular 
variation by downward continuation of the main field at the core mantle 
boundary (see for example Madden and Le Mouel, 1932). 
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TABLE 1 
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CHROFIN Corresponding date Number of processed days 

file number 


8A 

12/02-06 

8B 

06-10 

9A 

10-14 

9B 

14-18 

lOA 

18-22 


4 

4 

4 

4 

4 


lOB 

12/22-26 

llA 

26-30 

IIB 

30-01/03 

12A 

03-07 

12B 

07-11 


4 

4 

4 

4 

4 


13A 

13B 

14A 

14B 

15A 


01/11-15 
15-19 
— 19-23 
23-27 
27-31 


4 

4 

0.8 ^ 

4 

4 


15B 

16A 

16B 

17A 

17B 


01/31-02/04 

04-08 

08-12 

12-16 — 

16-20 


4 

4 

4 


0.8 4 

4 


18A 

18B 

19A 

19B 

20A 


02/20-24 

24-28 

— 28-03/03 

— 03-07 
07-11 


4 

4 

0.7 

0.2 ♦ 


4 


20B 

03/11-15 

21A 

15-19 

21B 

19-23 

22A 

23-27 

22B 

27-31 


0 « 

4 

4 

4 

4 
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TABLE 2 


Po 

av.RMSX 

av.RMSY 

av.RMSZ 

0 

23.5 

24.4 

17.7 

0.001 

23.6 

24.6 

17.8 

0.01 

24.2 

24.9 

19.2 

0.05 

24.9 

25.1 

19.9 

0.1 

25.3 

25.2 

19.9 


TABLE 3 


av.RMSX av.RMSY av.RMSZ 


MGSTOl 

20.6 

20.4 

14.2 

MGST02 

27.1 

30.0 

18.4 

MGST03 

20.0 

13.9 

13.6 


TABLE 4 : MGST03 coefficients 


1 

0 

-29961 


1 

1 

-1992 

5591 

2 

0 

-1998 


2 

1 

3040 

-2109 

2 

2 

1666 

-178 

3 

0 

1310 


3 

1 

-2178 

-348 

3 

2 

1247 

287 

3 

3 

838 

-236 

4 

0 

937 


4 

1 

781 

219 

4 

2 

401 

-251 

4 

3 

-419 

45 

4 

4 

206 

-291 

5 

0 

-186 


5 

! 

358 

49 

5 

2 

258 

152 

5 

3 

-72 

-152 

5 

4 

-160 

-82 

5 

5 

-51 

90 

6 

0 

50 


6 

1 

65 

-14 

6 

2 

41 

95 

6 

3 

-194 

66 

6 

4 

5 

-43 

6 

5 

13 

-2 

6 

6 

-108 

13 

7 

0 

99 


7 

1 

-58 

-81 

7 

2 

1 

-29 

7 

3 

20 

-5 

7 

4 

-13 

16 

7 

5 

0 

18 

7 

6 

12 

-23 

7 

7 

-1 

-10 

8 

0 

20 


8 

1 

6 

6 

8 

2 

1 

-18 

8 

3 

-11 

3 

8 

4 

-6 

-22 

8 

5 

4 

9 

8 

6 

2 

16 

8 

7 

7 

-13 

8 

8 

0 

-15 

9 

0 

25 


9 

1 

12 

-20 

9 

2 

0 

14 

9 

3 

-13 

8 

9 

4 

10 

-5 

9 

5 

-3 

-6 

9 

6 

-1 

9 

9. 

7 

6 

10 

9 

8 

1 

-6 

9 

9 

-5 

2 


TABLE 4 (continued) 



10 0 -1 

10 1 -5 

10 2 3 

10 3 -5 

10 4 -3 

10 5 5 

10 6 3 

10 7 ! 

10 8 2 

10 9 3 

10 10 0 

110 12 

111 0 

112 -2 

113 2 

11 < 0 

1’ 5 -1 

116 0 

117 2 

118 2 

11 9 -1 

IIP 2 

1 1 1 1 4 





0 

1 

1 

6 

-4 

0 

-1 

4 

0 

-6 

2 

0 

-2 

-3 

1 

0 

-2 
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